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Abstract —  New  architectures  and  an  associ¬ 
ated  CAD  method  are  proposed  for  linear  con¬ 
trollers.  The  architecture  is  suggested  by  re¬ 
cent  results  which  parametrize  all  controllers 
that  stabilize  a  given  plant. 

With  this  architecture,  the  design  ,  of  con¬ 
trollers  is  a  convex  programming  problem 
which  can  be  solved  numerically  even  for  large 
systems.  Constraints  on  the  closed-lobp  sys¬ 
tem  such  as  asymptotic  tracking,  decoupling, 
limits  on  peak  excursions  of  variables,  step  re¬ 
sponse  settling  time  and  overshoot,  as  well  as 
frequency  domain  inequalities  of  any  /fl  vari- 


conirol  or  actuator  inputs.  Its  outputs  we  de¬ 
compose  into  y,  the  measured  or  sensor  out¬ 
puts,  and  z,  the  regulated  variables.  Our  job 
is  to  design  a  controller  K,  with  input  y  only, 
and  output  u. 
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sponse  settling  time  and  overshoot,  as  Well  as 

frequency  domain  inequalities  of  any  //<J>  vari-  Controller 

ety  are  readily  incorporated  in  the  design.  The 

minimization  objective  is  general  and  can  in-  ^ 

elude  LQG,  //<*>,  new  t\  types,  or  any  combina-  ' 

tion  thereof.  Jl  .  d  _L 

To  specify  the  constraints  and  objectives  Figure  I:  Basic  plant  and  controller, 

mentioned  above,  we  are  developing  a  control 
specification  language.  *■  This  control  specification 

language  will  be  the  input  to  a  compiler  which  The  signal  y  represents  the  signal  actually 

will  translate  the  specifications  into  a  standard  accessible  to  the  controller  K,  including  any 

convex  program  in  It*-.  A  small  but  powerful  command  inputs,  which  may  be  considered  cx- 

subset  of  the  language  has  been  specified  and  ogenous  jnputs  (j.e.  components  of  tu)  passed 

its  associated  compiler  implemented.  The  re-  ,  r  .  r  2  c- 

„  .  directly  to  some  of  components  of  y/  Sim- 

sults  are  very  encouraging.  ,  • 

The  new  architecture  not  only  makes  design  lIarly  t}le  slSnal  u  presents  those  inputs  to 

of  the  controller  simple  but  also  its  tmplcmenia-  plant  which  our  controller  may  vary,  that 

tion.  These  controllers  can  be  built  right  now  is,  those  inputs  to  the  plant  manipulablc  by 


from  off  the  shelf  components  or  integrated  into 
silicon  using  standard  VLSI  cells. 

Details  will  appear  soon  in  a  full  length  pa¬ 
per. 

I.  Basic  Setup 

The  basic  plant  we  consider  is  shown  in  Fig¬ 
ure  1.  We  decompose  its  input  into  two  signal 
vectors,  tv,  the  exogenous  inputs,  and  u,  the 

*  Research  supported  in  part  by  ONR  under  NOOOM- 

8&-K-0112  and  NSF  under  ECS-85-52465. 


the  controller.  Thus  it  is  by  definition  that 
the  controller  has  input  y  only  and  output  u 
only. 

The  exogenous  input,  ui,  will  include 
real  physical  disturbances  or  noises  (torques, 
forces)  acting  on  the  plant,  any  actuator  or 
sensor  noise,  and,  as  mentioned  above,  any 

*We  may  think  of  the  exogenous  inputs  which  rep¬ 
resent  commands  as  the  actual  physical  commands 
inputs,  e.g.  angle  of  a  potentiometer,  and  the  corre¬ 
sponding  components  of  y  as  the  'sensed  command’, 
e.g.  voltage  output  from  a  potentiometer. 


command  inputs,  tv  may  also  contain  ficti¬ 
tious  inputs  injected  anywhere  in  the  plant. 

The  signal  z  represents  any  system  signal 
about  which  we  will  express  a  specification, 
regardless  of  its  accessibility  to  the  controller. 
Obvious  examples  are  the  actual  positions, 
forces,  temperatures,  etc.  we  wish  to  regulate 
or  control,  z  may  include  ‘internal  states’  or 
variables  we  wish  to  limit.  Some  components 
can  be  fictitious,  e.g.  linear  combinations  of 
state  variables  or  even  filtered  versions  of  sig¬ 
nals.  z  may  also  include  components  of  u  or 

y- 

Thus  Htw,  the  multi-input  multi-output 
closed-loop  map  from  w  to  z,  contains,  by  def¬ 
inition,  every  clostd-loop  map  of  interest. 

To  illustrate  this  decomposition  of  inputs 
and  outputs,  we  present  a  simple  example. 
A  motor  with  shaft  encoder  is  used  to  regu¬ 
late  the  angle,  0,  of  an  arm  to  some  (small) 
commanded  angle,  0rtj,  while  a  disturbance 
torque,  d,  acts  on  the  arm;  see  Figure  2.  Let 


(Motcr  Input  Vfcltage)  (Shaft  Encoder  Output) 

Fig  <i re  2;  Example  system. 

0,t  denote  the  shaft  encoder  output,  and  let 
n,e„»  =  0  —  0,e,  the  difference  between  the  ac¬ 
tual  angle  and  the  sensed  angle,  that  is,  sensor 
quantization  ‘noise’.  Exogenous  inputs  are 


There  are  two  sensed  outputs, 


Note  that  the  command  0rtj  is  simply  passed 
through  P;  thus  the  block  diagram  of  Figure 
1  could  be  drawn  in  a  more  conventional  way 
as  a  two  degree  of  freedom  SISO  controller. 

There  is  only  one  actuator  input,  the  motor 
input  voltage,  Vm;  that  is,  u  =  Vm.  A  possible 


choice  of  regulated  variables  is 


where  r  is  the  motor  torque.  We  will  refer  to 
this  example  several  times  in  the  sequel. 

We  assume  only  that  the  I/O  map  of  the 
plant,  P,  is  linear.  For  the  moment,  we  do 
not  say  whether  the  system  is  continuous  or 
discrete  time,  since  much  of  the  ensuing  dis¬ 
cussion  is  independent  of  this.  In  any  case,  a 
real  system  is  almost  certainly  hybrid  with  a 
continuous  time  physical  plant  and  a  discrete, 
possibly  multirate,  controller. 

We  partition  the  map  P  as 


Of  course  we  have 

H,w  =  Plw  +  PtuI< (/  -  Pyu/f )_1  Pyw  (1) 

where  H,w  is  the  closed-loop  map  from  the 
exogenous  inputs  w  to  the  regulated  variables 
z. 

Note  that  HIW  depends  on  the  controller 
K  in  a  linear  fractional  fashion.  A  relatively 
simple  constraint  on  Htw,  e.g.  that  a  certain 
entry  be  zero,  corresponds  via  (1)  to  a  much 
more  complicated  constraint  on  the  controller 
K.  We  remark  here  that  if  Pyo  =  0  then  (1) 
simplifies  considerably  to 

Htw  =  P tw  +  PtuKPyw- 

Here  II ,w  is  ajjine  in  K.  Note  that  this  cor¬ 
responds  to  the  case  where  there  is  essentially 
no  feedback  in  the  system. 

tj  Parameterization  of  Sta¬ 
bilizing  Controllers 

In  this  section  we  discuss  the  recent 
parametrization  of  all  stabilizing  controllers. 
Precise  definitions  and  all  details  can  be  found 
in  Vidyasagar’s  book  [1]. 

A  basic  requirement  is  that  the  controller 
stabilize  the  plant.  A  recent  major  theme  in 
control  theory  is  that  the  set  of  IItw’s  achiev¬ 
able  with  controllers  which  stabilize  the  sys¬ 
tem  is  a  linear  variety,  that  is,  a  translation  of 
a  linear  subspace. 


A  linear  variety,  C,  in  a  vector  space,  V,  is  I  =  the  identity  map.  A  little  calculation 
often  described  as  the  nullspacc  or  range  of  an  yields  the  Q-paramctrization  formula,  HiVJ  = 
afiine  map  from  or  to  V.  Thus  we  might  have  7\  +  T2QT3,  with 

C  —  {v  6  V  |  Av  —  6}  T\  —  Piw  A  P mDX  P yu), 


where  A  is  a  linear  map  from  V  to  W  and 
b  G  W.  This  is  a  description  of  C  in  terms  of  a 
linear  equality  constraint,  if  W  is  Rl ,  then  we 
can  interpret  each  entry  of  Av  =  b  as  a  single 
linear  functional  equality  constraint  in  V.  Let 
W  =  {Hzw  |  system  is  closed-loop  stable},  the 
set  of  achievable  Hzw's.  7i  may  be  described 
via  linear  equality  constraints  using  the  inter¬ 
polation  conditions  [6] 

Alternatively  we  may  describe  a  linear  vari¬ 
ety,  C,  as  the  range  of  an  affine  map: 

C  =  {Cm  +  d  |  u  €  U) 

where  C  is  a  linear  map  from  U  to  V  and  d  £  It- 
Such  a  representation  can  be  regarded  as  a 
free  parametric  representation  of  £;  u  is  a  free 
parameter. 

A  free  parametric  representation  of  7i  can 
be  derived  using  stable  coprime  fractional 
theory  [6,1]  One  standard  form  is  the  ‘ Q - 
parametrization’: 

n  =  {T!+  T2QT3  I  Q  stable}. 

where  T\,  T2,  and  T3  are  some  stable  maps. 

Note  that  the  Q-paramctrization 

ff,u,=Tl+T2QT3  (2) 

has  exactly  the  same  form  as  the  no-feedback 
formula  mentioned  in  Section  1;  we  will  soon 
see  a  block  diagram  interpretation  of  the  Q- 
parametrization  in  which  the  parameter  Q 
‘sees’  no  feedback.  In  particular,  Htw  is  affine 
in  Q. 

The  derivation  of  the  Q-parametrization 
starts  with  any  controller,  Ka,  which  stabilizes 
the  plant  and  which  we  will  call  the  nominal 
controller.  Let  I\0  =  be  a  left  stable 

coprimc  factorization  [1]  of  the  nominal  con¬ 
troller  and  Pyu  =  D~lN  be  a  left  stable  co¬ 
prime  factorization  of  Pyu.  Then  every  con¬ 
troller  of  the  form 

I<  =  {Y  -  QN)~l(X  +  QD),  Q  stable 

stabilizes  P  and  conversely  every  controller 
which  stabilizes  P  has  this  form  for  some  sta¬ 
ble  Q.  There  is  a  similar  characterization  of 
the  stabilizing  controllers  in  terms  of  right  co- 
prime  factorizations. 

Since  A'o  stabilizes  Pyu,  there  is  a  right  co¬ 
prime  factorization  of  Pyu-  with  XN  +  YD  = 


T2  =  PiuD, 

T3  =  DPyW. 

A  simple  block  diagram  of  the  Q- 
parametrization  is  shown  in  Figures  3  and  4. 


Figure  3:  Block  diagram  of  Q-paramctrization. 


Figure  4:  Q-paramctrization  as  modification  to  nominal 
controller. 

7}  is  simply  the  Htw  achieved  with  the  nom¬ 
inal  controller  I\0\  T2  is  the  map  from  w  to  e, 
T3  is  the  map  from  v  to:  (see  Figure  4).  The 
key  to  the  parametrization  is  that  the  closed- 
loop  map  from  v  to  c  is  zero,  so  that  Q  ‘secs’ 
no  feedback. 

Doyle  [7]  has  given  a  very  nice  interpreta¬ 
tion  of  the  Q-paramctrization  when  the  nom¬ 
inal  controller  is  an  estimated  state  feedback. 
In  this  case  e  is  simply  the  output  prediction 
error,  y  —  y,  of  the  observer,  and  t;  is  just 
added  before  the  observer  tap  to  the  output  of 
the  nominal  controller  as  shown  in  Figure  5. 
The  controller  shown  in  Figure  5  is  sometimes 


called  an  observer-based  controller  (ODC);  the 
point  is  that  every  controller  which  stabilizes 
P  can  be  realized  as  an  observer-based  con¬ 
troller. 


Figure  5:  Doyle's  interpretation  of  Q-parametrization 
for  estimated  state  feddbeak  nominal  controller.  Here 
F  is  a  stabilizing  state  feedback  gain  and  y  =  C x  where 
x  is  the  plant  state. 

m  Closed-Loop  Constraints 
and  Objectives 

In  this  section  and  the  next,  we  consider  the 
problem  of  designing  the  parameter  Q,  given 
the  nominal  maps  Tlt  T2,  and  T3.  We  first  ob¬ 
serve  that  many  typical  requirements  on  the 
closed-loop  performance  of  our  system  result 
in  convex  constraints  on  II  xw  and  thus  on  the 
parameter  Q.  Before  proceeding,  we  note  that 
requirements  on  the  open-loop  system  gener¬ 
ally  do  not  result  in  convex  constraints  on  IIXW 
or  Q.  One  important  example  is  the  require¬ 
ment  that  the  controller,  K  be  diagonal,  that 
is,  that  I\  be  a  decentralized  controller.  An¬ 
other  important  example  is  the  requirement 
that  the  controller  be  stable. 

Suppose  our  design  problem  is  specified  as 

minimize  $(IIXW) 

//hi  6  A. 

cl.  loop  stable 

where  <I>  is  a  convex  functional  and  K,  a  convex 
(constraint)  set.  This  is  equivalent  to 

minimize  ^'(Q) 

Q  g  )C  W' 

where 

ic  =  {Q  |  Ti  +  T2QT3  e  c} 

and 

4>(Q)  =  <I>(7\  +  T2QT3). 


Note  that  <I»  and  K.  are  convex.  Of  course 
K,  may  be  empty;  this  simply  means  that  no 
stabilizing  controller  can  satisfy  the  constraint 
H zxu  G  L. 

We  now  list  some  typical  constraints  (and 
objectives)  on  Hxw,  some  collection  (sum)  of 
which  might  describe  the  constraint  set  K.  (ob¬ 
jective  functional  $). 

For  purposes  of  discussion,  we  assume  the 
system  is  discrete-time.  H  will  denote  the 
transfer  matrix  of  IIXW,  h  its  impulse  response 
matrix,  and  s(t)  =  5Z<=o  MO  *ts  steP  re" 
sponse.  Of  course  a  convex  constraint  or  func¬ 
tional  of  H ,  h,  or  s  is  a  convex  constraint  or 
functional  of  HIW. 

Where  possible,  we  will  refer  to  the  example 
system  of  §1,  which  we  presume  has  been  sam¬ 
pled  at  some  appropriate  rate.  When  referring 
to  this  system,  we  will  use  symbolic  subscripts 
written  in  square  brackets  so  the  reader  need 
not  refer  to  §1  to  find  which  components  of 
tv  and  z  have  which  interpretations.  Thus  we 
will  write  H[0][0TeJ}  instead  ol  IIu. 

•  Asymptotic  Tracking,  Decoupling,  and 
Regulation  - 

The  step  response  from  some  command 
input  to  the  regulated  variable  it  is  sup¬ 
posed  to  control  must  converge  to  one, 
e.g- 

sWr«/K0  =  L 

Equivalently, 

JWr./)(e>°)  =  1. 

This  constraint  is  a  single  linear  func¬ 
tional  equality  constraint  on  II,  hence, 
of  course,  a  convex  constraint. 

Asymptotic  tracking  of  ramps  or  more 
complicated  inputs  can  be  handled  as 
two  or  more  linear  functional  equality 
constraints,  e.g.  //[0l[0re/](M°)  =  1, 

ir[o][dr'j](c>°)  =  0. 

Asymptotic  regulation  and  asympolotic 
decoupling  arc  similar  constraints.  We 
may  require  that  a  regulated  variable 
asymptotically  reject  constant  inputs 
appearing  at  certain  cxogcncous  in¬ 
puts.  When  the  exogenous  input  is  an¬ 
other  command  input,  this  constraint  is 
asymptotic  decoupling;  when  the  exoge¬ 
nous  input  is  some  disturbance,  this  con¬ 
straint  is  asymptotic  regulation.  For  ex¬ 
ample,  to  ensure  that  0  asymptotically 


rejects  any  constant  disturbance  torque, 
we  might  specify 

H[0  -  Orcj)[d}(ei°)  =  0. 

•  Overshoot,  Undershoot,  and  Settling- 
Time  Limits 

We  may  require  that  some  step  response 
lie  in  the  unshaded  region  of  Figure  6, 

eg- 

-0.3  <  s[0][0re/](O  <  1-3;  0  <  t  <  10, 

l«MPr./]- i(0l  <V(<  +  i);  <>io- 

This  constraint  can  be  expressed  as  a 
collection  of  linear  functional  ineqqquali- 
ties  on  s  (hence  II 2W)'  L(t )  <  S32  <  U(t) 
for  t  =  0, 1, ....  Thus  the  set  of  II2W ’s  sat* 

•  isfying  this  constraint  is  convex  (it  may 
of  course  be  empty). 


Figure  6:  Step  response  overshoot,  undershoot,  and  set¬ 
tling  time  constraints. 

•  Bounds  on  Closed-Loop  Signal  Peaks 
Given  bounds,  IV j,  on  each  exogenous 
input,  we  may  require  that  each  regu¬ 
lated  variable  be  bounded  by  some  given 
maximum  Z3.  This  constraint  could 
arise  from  the  requirement  not  to  sat¬ 
urate  an  actuator  or  sensor  or  exceed 
some  internal  variable  force,  torque,  or 
current  limit.  This  constraint  is  equiva¬ 
lent  to 

^»»B|  OO 

Y  Wi  Yi  M0I  <  for  *  =  1.  -Nr eg 
;  =  1  1=0 

where  Ntxag  is  the  number  of  exogenous 
inputs  and  Nreg  the  number  of  regulated 
variables. 

Thus  the  constraint 

OO 

1=0 

+  o.osf;iMi/«MOI 

1=0 

co 

+  0.001  £>[V’mM0|<  50 

1  =  0 


will  guarantee  that  whenever  the  com¬ 
mand  input  0ref  is  bounded  by  0.1,  the 
disturbance  d  by  0.03,  and  the  sensor 
noise  n,en3  by  0.001,  the  motor  voltage 
input  Vm  will  be  bounded  by  50. 

Again,  this  is  a  convex  constraint  on  h, 
and  thus  on  H 2W. 

•  Small  RMS  Disturbance  Response 
If  the  exogenous  input,  w,  is  driven  by 
a  wide-sense  stationary  stochastic  pro¬ 
cess  with  some  specified  spectral  distri¬ 
bution,  then  E :(t)TGTG:(t),  where  G 
is  some  weighting  matrix,  should  be  as 
small  as  possible.  Typically,  the  stochas¬ 
tic  process  would  have  zero  power  in 
the  exogeneous  input  channels  which  are 
command  inputs,  and  G  would  ha%-e  only 
a  few  nonzero  entries. 

We  have 

E z(t)TGTGz(t)  =  x  the  Trace  of: 

G  J*'  H(ein)Su,(e*n)H(e~in)TdQ  GT 

where  Sw  is  the  power  spectral  den¬ 
sity  of  the  stochastic  process  driving  w. 
This  objective  is  a  nonnegative  quadratic 
functional  of  II  and  thus  of  H2W ,  in  par¬ 
ticular  it  is  convex.3 
For  our  example,  a  possible  Sw  would  be 

'  0  0  0  1  (0rt/) 

Sw(dn)=  0  S„(e;n)  0  (d) 

0  0  A:/12 

where  5j(e;n)  is  the  spectral  density  of 
the  disturbance  torque  and  A  is  the  step 
size  of  the  shaft  encoder.  Note  that  we 
put  no  noise  power  in  the  command  in¬ 
put  0rej.  If  we  arc  interested  in  the  noise 
power  in  the  pointer  angle,  0,  then  we  let 
G  =  [  1  0  0  0  ]  and  the  functional  above 
is  simply 


\H[0M(cin)\iSd(cin)dQ 


3llcre  quadratic  functional  includes  a  linear  function 
and  a  constant  (functional)  term. 


which  can  be  interpreted  as  the  total 
noise  power  in  0  due  to  the  disturbance 
and  sensor  noise. 

Bounds  on  Transfer  Function  Peak  Mag¬ 
nitudes 

We  may  specify  an  upper  bound  for  some 
entry  (or  block  of  entries)  of  //,  e.g. 
|tf6s(ein)|  <  U(Q)  for  all  Q,  where  U(Q) 
is  some  bound  function.  This  constraint 
may  arise  in  several  ways. 

First,  many  classical  specifications  of 
bandwidth  or  peaking  are  expressed  this 
way.  Referring  to  our  example,  we  might 
specify 

\H[0-0rtJ)[0r'j](ein)\  <  0.03  for  |fl|  <  VB 

so  that  for  frequencies  below  CtB,  we  get 
-30db  of  rejection. 

\H[0  -  0M/)[0r./](ein)|  <  3  for  all  fi 

to  keep  gain  in  the  command  to  track¬ 
ing  error  map  under  about  lOdb  for  all 
frequencies. 

The  second  way  a  transfer  function  peak 
magnitude  constraint  can  arise  is  as  fol¬ 
lows.  If  we  require  a  small  RMS  response 
of  some  regulated  variable,  as  in  the  pre¬ 
vious  constraint,  but  know  only  a  bound 
on  the  total  power  in  the  process  driv¬ 
ing  the  exogenous  input  w,  and  not  its 
spectral  distribution,  then  we  are  led  to 
a  transfer  function  peak  magnitude  con¬ 
straint;  this  is  the  basis  of  IIoo  control 
theory. 


Such  a  constraint  can  be  used,  for  ex¬ 
ample,  to  guarantee  stability  of  the  sys¬ 
tem  with  a  saturating  actuator.  All  con¬ 
straints  involving  step  responses  can  be 
generalized  to  other,  arbitrary,  input  sig¬ 
nals,  or  even  sets  of  signals.  Thus  we 
may  specify  a  maximum  peak  tracking 
error  for  command  inputs  bounded  by 
Bcmd  and  slew  rate  under  RCmi- 

Finally  we  mention  that  the  boundary  be¬ 
tween  constraint  and  objective  is  not  sharp;  we 
could,  for  example,  try  to  minimize  an  over¬ 
shoot  in  a  certain  step  response,  or  put  hard 
constraints  on  the  RMS  response. 

At  this  point  it  should  be  clear  that  a  large 
number  of  practical  constraints  on,  and  objec¬ 
tives  for,  closed-loop  system  performance  can 
be  formulated  as  a  standard  convex  program 
for  the  parameter  Q. 

This  program  for  Q  is  infinite  dimensional 
and  generally  cannot  be  solved  analytically  ex¬ 
cept  in  special  cases.  The  small  RMS  distur¬ 
bance  response  objective  alone  can  be  solved 
using  Wiener-IIopf  or  LQG  methods;  certain 
forms  of  the  small  peak  transfer  function  ob¬ 
jective  problem  are  solved  by  Iloo-optimal  con¬ 
trol  theory,  certain  forms  of  the  small  peak  dis¬ 
turbance  objective  problem  are  solved  by  the 
new  ^-optimal  control  theory  introduced  by 
Vidyasagar  [S]  and  developed  by  Dahlch  and 
Pearson  [2]. 


IV. 


Specifications  to  Convex 
Program  Compiler 


Third,  a  sufficient  condition  to  stabilize 
not  just  the  plant,  but  all  plants  ‘near’ 
our  given  plant  (in  a  peak  frequency  re¬ 
sponse  deviation  sense),  may  be  formu¬ 
lated  as  a  bound  on  certain  closed-loop 
frequency  response  magnitudes. 

These  IIoo  type  constraints  are  convex 
constraints  on  II  and  hence 


The  translation  from  closed-loop  constraints 
or  objective  on  H,w  to  convex  program  for 
the  parameter  Q  is  mathematically  straight¬ 
forward.  As  mentioned  above,  the  resulting 
program  is  infinite  dimensional.  An  approxi¬ 
mate  solution  can  be  found  numerically  when 
the  parameter  Q  is  restricted  to  a  large,  but 
finite,  dimensional  space.  Let 


Miscellaneous  Bounds  and  Objectives 

We  mention  here  some  less  common  con¬ 
straints.  A  slew-rate  constraint  on  the 
step  response,  say  sn,  may  be  enforced 
as: 

|/»n(OI  <  R  for  all  t. 

A  passivity  or  minimum  dissipation  con¬ 
straint,  say  on  the  1,  1  entry  of  IIIW,  can 
be  formulated  as 

3i//u(cin)  >  -D  for  all  0. 


Q  =  £*.r. 

1=1 

where  x.  €  R  and  T,  arc  fixed  stable  maps. 
We  will  call  x  -  [xi , Xl]t  the  decision 
variables.  With  this  additional  restriction,  we 
have  a  standard  finite  dimensional  convex  pro¬ 
gram  in  Rl 

minimize  /(x)  ... 

xe)C  '  '  ' 


While  the  convex  constraint  set  AC  cannot  be 
explicitly  found,  it  is  easily  approximated,  and 
then  this  explicit  finite  dimensional  convex 
program  can  be  solved  numerically  by  any  of 
a  number  of  methods. 

Let  us  first  briefly  discuss  the  choice  of 
parametrization  and  the  r,’s,  since  this  aflccts 
how  well  the  finite  dimensional  program  (4) 
approximates  the  program  (3). 

In  theory,  it  makes  no  difference  which 
nominal  controller  nor  which  coprime  factor¬ 
ization  is  used  to  form  the  Tfs  in  the  Q- 
parametrization;  all  parametrizations  produce 
precisely  the  same  set,  7i ,  of  Hzw' s.  The  par¬ 
ticular  Q  which  yields  a  given  Hlw,  however, 
will  depend  on  the  choice  of  Tfs.  As  men¬ 
tioned  above,  the  Q’s  will  be  restricted  to 
some  finite  dimensional  space,  so  the  set  of 
Ht u,’s  achievable  in  practice  is  a  finite  dimen¬ 
sional  linear  variety  containing  T\,  which  is 
the  Hzw  achieved  with  the  nominal  controller. 
This  means  that  the  nominal  controller  used 
should  yield  at  least  reasonable  performance, 
so  that  we  have  more  of  a  chance  of  finding 
a  ‘moderate’  Q  which  achieves  good  perfor¬ 
mance.  By  ‘moderate’  we  mean  a  Q  which 
is  neither  too  ‘big’  nor  too  ‘stiff,  that  is,  has 
too  broad  a  spread  of  dynamics.  An  inappro¬ 
priate  nominal  stabilizing  controller  will  place 
severe  demands  on  the  filter  Q.  For  exam¬ 
ple,  a  state  feedback  which  places  the  eigen¬ 
values  far  to  the  left  of  the  bandwidth  achiev¬ 
able  when  other  factors  such  as  control  au¬ 
thority  are  taken  into  account,  will  require  a 
‘large’  Q  filter,  essentially  to  cancel  the  unrea¬ 
sonably  large  plant  input  signal  generated  by 
the  excessive  state  feedback.  So  in  practice, 
the  choice  of  nominal  controller  is  not  arbi¬ 
trary. 

A  similar  comment  applies  to  the  choice  of 
coprimc  factorizations,  which  will  affect  the 
maps  Tn  and  T3.  Roughly  speaking,  they 
should  be  of  reasonable  size  and  have  most 
of  their  dynamics  not  too  far  from  the  final 
closed-loop  bandwidth. 

We  turn  now  to  the  practical  prob¬ 
lem  of  translating  a  large  number  of  con¬ 
straints/objective  on  step  responses,  transfer 
functions,  ll.MS  responses,  and  so  on,  into  the 
corresponding  convex  program  for  the  decision 
variables.  We  propose  that  this  job  be  mech¬ 
anized  by  the  use  of  a  compiler  which  accepts 
as  input  a  list  of  closed-loop  constraints  and 
a  description  of  the  objective  and  as  output 


produces  the  convex  program  for  the  decision 
variables.  The  constraints  and  objective  would 
be  specified  in  a  control  specification  language, 
which  would  be  natural  for  the  control  engi¬ 
neer,  refering  directly  to  step  responses,  noise 
powers,  transfer  functions,  etc.  The  output  of 
the  compiler  would  be  a  convex  program  which 
could  be  solved  by  a  convex  program  solver,  as 
shown  in  Figure  7. 

I - T - 1 


Figure  7;  Design  of  appro xiniatcly  optimal  controller 
using  control  specification  to  convex  program  compiler. 

IV  A  Control  Specification 
Language 

Here  we  discuss  the  control  specification,  the 
input  to  the  compiler.  Many  of  the  features  we 
discuss  have  actually  been  implemented  in  our 
compiler;  §Vl  describes  what  wc  have  imple¬ 
mented  so  far. 

An  overall  general  structure  for  the  control 
specification  is 

constraints  { 

<constraint>; 

<constraint>; 


« 


objective  { 

description  of  objective 

> 

For  example,  the  constraints  might  have  the 
forms  encountered  in  §111,  such  as 
constraints  { 

=  1; 

-0.3  <  s{0][0rtJ){i)  <  1.3  for  0  <  t  <  10; 

0-lE“o  IMV'mJPr./KOI 
+0.03  £~0  IMi'mMOl 
+0.001  2“olMV’mM0l<  50; 

3t/Iu(ejU)  >-D  for  all  fi; 

} 

A  description  of  the  objective  could  be 
objective  { 

&/o'raM(^n)laWn)^ 

For  a  detailed  description  of  the  syntax,  see 
§VI. 

v  Two  Proposed  Controller 
Architectures 

Let  Q‘  =  £f_j  *‘iQi  be  a  satisfactory  Q  for 
the  control  problem.  The  resulting  controller 
A'*  =(Y  -Q'N)~l(X  +  Q'D) 
will  generally  be  full,  (each  input  will  alTcct 
every  output)  and  have  a  large  number  of 
states,  perhaps  hundreds.  The  standard  pro¬ 
cedure  at  this  point  would  be  to  apply  some 
sort  of  model  reduction  to  arrive  at  a  low  or¬ 
der  controller,  Kred,  with  only  tens  of  states. 
We  propose  instead  two  methods  for  imple¬ 
menting  the  full  A*,  both  based  on  the  Q- 
parametrization.  Doth  methods  rely  on  the 
availability  and  ease  of  VLSI  implementation 
of  many-tap  FIR  filters,  used  extensively  in 
signal  processing,  for  example  in  equalizers  for 
modems.  Thus,  our  controllers  will  tend  to 
have  far  greater  orders  than  is  currently  com¬ 
mon.  While  it  is  certain  that  there  is  a  much 
lower  order  I\red  which  provides  satisfactory 
performance,  we  question  the  need  to  find  it. 

Modified  Nominal 
V.A  Controller  Architec¬ 
ture 

Consider  the  block  diagram  of  the  Q- 
parametrization  given  in  Figure  4.  We  pro¬ 
pose  that  controllers  actually  be  built  with  the 


architecture  shown  in  Figure  4,  with  Q  real¬ 
ized  (possibly)  separately  in  hardware  as  an 
FIR  filter.  This  architecture  can  be  thought  of 
as  a  simple  modification  of  the  nominal  con¬ 
troller;  the  nominal  controller  must  be  modi¬ 
fied  to  yield  an  auxiliary  output  signal  e  and 
accept  an  auxiliary  input  signal  v.  In  some 
controllers,  the  signal  e  is  easily  available,  even 
though  Figure  3  suggests  that  we  need  an  addi¬ 
tional  N  and  D.  Similarly,  injecting  the  signal 
v  as  shown  in  Figure  4  is  also  straightforward. 

For  example,  if  the  nominal  controller  is  an 
estimated  state  feedback,  then  as  mentioned 
earlier,  e  is  the  output  prediction  error  of  the 
observer,  and  v  is  just  added  to  the  output 
of  the  nominal  controller,  before  the  observer 
tap,  as  shown  in  Figure  5.  Thus,  in  many  cases 
the  hardware  and/or  computational  costs  of 
modifying  the  nominal  controller  are  slight, 
since  the  Q  filter  does  all  of  the  additional  pro¬ 
cessing. 

We  note  that  this  method  can  be  used  to 
improve  existing  controllers. 

v  B  FIR  Coprime  Archi¬ 
tecture 

Generally,  the  nominal  controllers  which 
provide  easy  access  to  auxiliary  output  e  are 
those  with  complexity  on  the  order  of  the 
plant,  the  estimated  state  feedback  nominal 
controller  a  good  example.  The  reason  is, 
roughly  speaking,  that  e  is  the  pseudo-state 
observer  error,  which  means  that  the  modified 
nominal  controller  really  contains  a  pseudo- 
state  observer  within  it,  which  is  likely  to  have 
the  same  complexity  as  the  plant.  This  is  ad¬ 
mittedly  vague,  and  only  ofTcrcd  as  our  inter¬ 
pretation. 

In  nominal  controllers  of  much  lower  com¬ 
plexity  than  the  plant,  e.g.,  a  diagonal  PI  con¬ 
troller,  e  is  usually  not  easily  accessible,  and 
must  be  reconstructed.  In  this  case  we  propose 
that  the  entire  controller  be  implemented  us¬ 
ing  FIR  filters. 

Consider  the  coprime  factorization  of  the 
controller  A'*  given  above.  Suppose  Y ,  Q’ , 
N,  A',  and  D  are  all  FIR,  or  very  nearly  ap¬ 
proximated  by  FIR  filters.  Recall  that  unlike 
the  plant  or  controller,  these  operators  are  sta¬ 
ble  and  hence  their  impulse  matrices  decay  to 
zero.  With  proper  choice  of  sample  rate  and 
parametrization,  100  tap  FIR  fillers  should  be 
more  than  adequate.  In  this  case,  A'*  can  be 


realized  as  a  feedback  connection  around  one 
FIR  filter,  as  we  now  show.  A'*  is  governed  by 
the  equations 

Yu  —  Xy  =  v  =  QXu  +  QD 
(see  Figure  3).  We  rewrite  these  equations  as 
ti  =  (I  -  V)u  +  Xy  +  v 


v  =  QXu  +  QDy 

which  we  interpret  as  an  FIR  filter,  F,  with 
input  u,  y  and  v,  and  output  u  and  v,  and 

F-\  X-  M 

QN  QD  0 

This  is  shown  in  Figure  8.  Of  course  there  arc 


Figure  6:  FIR  coprimc  iucliilccturc  for  controllers. 


many  other  F’s  which  realize  A'*,  and  we  do 
not  vet  know  a  good  procedure  for  picking  one. 

Finally,  we  mention  that  for  both  archi¬ 
tectures,  it  is  possible  to  vary  the  effect  of 
the  parameter  Q  from  ‘off  to  ‘on’  by  insert¬ 
ing  a  gain,  o,  which  varies  between  zero  and 
one,  in  the  right  place  (e.g.  in  the  bottom 
loop  of  Figure  8).  Thus  we  have  a  homo - 
iopy  (continuous  deformation)  of  our  nomi¬ 
nal  controller  A'o  into  our  designed  controller 
K' ,  with  all  intermediate  controllers  stabi¬ 
lizing  the  plant,  indeed  with  Hlw  affine  in 
the  parameter  a.  Such  a  parameter  might 
be  quite  useful  as  a  control  authority  param¬ 
eter.  If  our  nominal  controller  has  somewhat 
lower  authority/bandwidth  than  our  designed 
controller,  but  is  more  robust  to  plant  vari¬ 
ations  or  failures,  then  a  can  vary  the  con¬ 
trol  authority  and  robustness  of  our  controller 
from  nominal  (resp.,  low  and  high)  to  designed 
(resp.,  high  and  low).  In  an  emergency,  we  can 
smoothly  drop  back  to  the  nominal  controller. 
Of  course,  this  will  only  work  with  moderate 
plant  variations,  and  cannot  take  the  place  of 
a  real  fault  tolerant  control  system  (with  fault 
detection  and  reconfiguration). 


VI. 


Implementation  of  a  Sun- 


SET 


We  have  implemented  a  compiler  which  rec¬ 
ognizes  a  small  but  powerful  subset  of  control 


specification  language.  Time  and  frequency 
domain  linear  and  magnitude  frequency  in¬ 
equality  constraints,  and  quadratic  objective 
arc  recognized.  The  system  must  have  only 
one  actuator,  and  the  filter  Q  is  simply  a  FIR 
filter.  We  are  currently  implementing  a  much 
more  general  compiler. 

We  assume  Q  has  n1?  inputs  and  a  single 
output-,  so  far  we  have  implemetcd  a  compiler 
for  control  problems  that  have  only  one  actu¬ 
ator  or  control  input.  Q  is  assumed  to  be  an 
iW/ytap  FIR  filter  in  every  channel.  Thus, 
there  are  a  total  of  n,,  x  nCOejj  parameters 
in  the  filter  Q  and  these  are  the  decision  vari¬ 
ables.  In  this  case 

n., 

Hlw  =  Tt  +  ^Q.Gi 

1=1 

where  Q,  is  the  map  from  the  ith  input  of  Q 
to  its  output  and 

Gi  =  T,cjT3 


(e,-  is  the  ith  unit  vector).  Since  Q  is  FIR 
with  ncoejj  taps  in  each  channel,  the  system 
impulse  response  at  a  given  discrete  time  k  is 

^ 1  4  Kcon* 

hiw{k)  =  t!(k)  +  J2  ^  <l,{j)g,{k  -  j), 

»=i  j  =o 

where 

riconv  —  min(ncoeyy  \,kf 

Clearly  the  system  impulse  response  at  any 
discrete  instant,  k,  is  affine  in  the  decision  vari¬ 
ables,  <7,0'). 

Our  compiler  recognizes  three  time  domain 
functionals: 

uste;<[i][j](F), 

M'lUK*). 

?[«](*)• 

These  arc,  respectively,  the  step  and  impulse 
responses  from  input  j  to  output  i  at  time  k, 
and  the  step  response  of  Q  from  it’s  input  i  at 
time  k  (recall  that  Q  has  a  single  output). 

Three  frequency  domain  functionals  arc  rec¬ 
ognized: 

AeJ/tOLilM), 

ImJl[i][j}(r,0), 

MagJl[i]{j}{r,0). 

These  arc,  respectively,  the  real  part,  imagi¬ 
nary  part,  and  magnitude  of  the  frequency  re¬ 
sponse  from  input  j  to  output  i  at  rc1$ ,  r  and 


0  are  real.  Since  our  current  convex  program¬ 
ming  package  allows  only  linear  constraints, 
constraints  involving  Magll  arc  parsed  into 
two  contraints  on  the  real  and  imaginary  parts 
of  // —  at  the  expense  of  at  most  a  l.odD  er¬ 
ror.  It  could  be  modified  to  parse  into  more 
constraints,  yielding  less  error. 

Real  scalar  variables  and  scalar  expressions 
involving  standard  mathematical  functions  are 
allowed. 

Functional  inequalities  are  allowed  to  take 
the  form 

/  <  e; 

/  >  e; 

/  =  e; 

c  <  /  <  e; 

I  /  I  <  e; 

where  /  represents  one  of  the  above  function¬ 
als  and  e  represents  any  scalar  expression. 

A  quadratic  objective  can  be  specified  as  the 
sum  of  any  number  of  the  following  terms: 

h.sqr{>\{]}{k), 

MagJljsqr[i\[j](r,0), 

<7-*7  »•['](*)• 

These  arc,  respectively,  the  square  of  the  im¬ 
pulse  response  from  input  j  to  output  i  at  time 
k,  the  square  of  the  frequency  response  mag¬ 
nitude  from  input  j  to  output  i  at  re*9 ,  and 
the  square  of  the  impulse  response  of  Q  from 
its  input  i  to  it’s  output  at  time  k. 

Nested  looping  is  allowed.  This  makes 
repetitive  constraints  and  objective  terms  easy 
to  specify,  for  example 

cor-traints  { 
for  t  =  0  to  10 

-  0.3  <  usfc;.[0)[0rt/](()  <  1.3; 

} 

objective  { 
for  fl  =  0  to  t  step  0.03 
0.03  .  ,U<ii/_//_s</7-[d][</j(c>n); 

} 

The  standard  C  language  preprocessor  c.pp 
is  used,  which  allows  the  use  of  macros  and  file 
include  facilities. 
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